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ABSTRACT 


The study focuses on the axisymmetric deformation response of unsymmetrically lami- 
nate cylinders loaded in axial compression by known loads. A geometrically nonlinear anal- 
ysis is used. Though buckling is not studied, the deformations can be considered to be the 
prebuckling response. Attention is directed at three 16 layer laminates: a (90,/0„) T ; a 

(0»/90,) t ; and a (0/90) a . The symmetric laminate is used as a basis for comparison, while the 
two unsymmetric laminates were chosen because they have equal but opposite bending- 
stretching effects. Particular attention is given to the influence of the thermally-induced pre- 
loading deformations that accompany the cool-down of any unsymmetric laminate from the 
consolidation temperature. Simple support and clamped boundary conditions are considered. 
It is concluded that: (1) The radial deformations of an unsymmetric laminate are significantly 
larger than the radial deformations of a symmetric laminate. For both symmetric and unsym- 
metric laminates the large deformations are confined to a boundary layer near the ends of the 
cylinder; (2) For this nonlinear problem the length of the boundary layer is a function of the 
applied load; (3) The sign of the radial deformations near the supported end of the cylinder 
depends strongly on the sense (sign) of the laminate asymmetry; (4) For unsymmetric lami- 
nates, ignoring the thermally-induced preloading deformations that accompany cool-down re- 
sults in load-induced deformations that are under predicted; and (5) The support conditions 
strongly influence the response but the influence of the sense of asymmetry and the influence 
of the thermally-induced preloading deformations are independent of the support conditions. 


INTRODUCTION 

Composite cylindrical structures are known for their efficiency. In addition, they are well 
suited for fabrication by automated fiber-placement techniques such as filament winding. 
Axial loadings are common, as are pressure, and to a lesser extent, torsional and bending 
loads. Cylinders subjected to axial end loads are susceptible to buckling or collapse if the 
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load reaches significant levels, or if the cylinder contains imperfections in the form of out-of- 
roundness, variations in thickness, material flaws, or combinations of these. In addition to 
their structural efficiency, composite cylindrical structures are somewhat more tolerant of un- 
symmetric lamination sequences than flat laminates. Despite the tendency of an unsymmetric 
laminate to deform significantly when cooled from its processing temperature, cylindrical 
forms must remain cylindrical due to the axisymmetric nature of the basic geometry. In many 
applications involving cylinders with hundreds of layers in the cylinder walls, the stacking ar- 
rangement is often unsymmet'ical simply because it is easier to continue a winding sequence 
through the wall rather than reverse the sequence half way through the wall. For the many- 
layered case, the degree of asymmetry is small and is not generally accounted for in design 
and analysis procedures. 

There is the general feeling that having any degree of asymmetry is not good. However, 
there is an increasing interest in unsymmetric laminates for structural applications. Unsym- 
metric laminates may indeed be the minimum weight design, the addition of more layers to 
make the laminate symmetric simply adding weight to the structure with no increase in per- 
formance. Also, manufacturing issues, just mentioned, may dictate that an unsymmetric 
laminate is less expensive to fabricate. And finally, there is a considerable degree of elastic 
coupling in unsymmetric laminates that is not present in symmetric laminates. For structural 
tailoring, these couplings may prove beneficial. This reports details the finding of a very basic 
study focused on determining some of the characteristics of unsymmetrically laminated cyl- 
inders, in particular, unsymmetrically laminated cylinders subjected to axial compressive 
loads. Considerable work has been done by others in the area of axially-loaded composite 
cylinders. However, in numerical examples in the past work that have involved unsymmetric 
laminates, no attention has been given to the fact that when an unsymmetrically laminated 
cylinder is actually constructed, the temperature change from the processing temperature to, 
say, room temperature results in a cylinder that has a deformed shape before it is even loaded 
[1-3]. That is, there is a preloading deformed shape that, when loaded, deforms further. When 
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considering prebuckling deformations in a buckling analysis, for example, or the response in 
a collapse analysis, does inclusion of these thermally-induced preloading deformations have 
an effect on the predicted result? The purpose of the work here is to begin to accurately 
model this thermally-induced deformation in unsymmetric laminates and determine if the 
thermally-induced preloading deformations, in addition to the presence of the bending- 
stretching of unsymmetric laminates, 1) has an influence on response, and; 2) can be used to 
advantage. As an example of point 2, it might be possible that bending-stretching coupling 
could be used to increase the buckling load, or alter the sensitivity of the cylinder to post- 
buckling collapse. 

The particular work reported on here focuses on a derivation of the equations governing 
the prebuckling deformations, including thermal effects, and the equations governing buckling. 
For the present, only perfect cylinders are considered, i.e., perfectly round, uniform material 
properties, etc. Buckling and collapse are not addressed. The next sections define the ge- 
ometry, coordinate system, and nomenclature used to study cylinder response, and then pro- 
ceed to derive the equilibrium and buckling equations. Though the latter are not used here, 
they are derived for future use. The derivation is based on variational principles and total 
potential energy. The first variation of the total potential energy is used to establish the 
equilibrium conditions, while the second variation, through the Trefftz stability criterion, is 
used to develop the buckling equations. Following the derivation the axisymmetric 
prebuckling responses of several simple unsymmetrically laminated cylinders to an axial 
compressive load are studied. The influence of the thermally-induced preloading defor- 
mations, and the influence of changing the sign of the bending-stretching coupling effects are 
studied. This prebuckling study is limited in that only three different cylinders are studied. 
These are: a (0,/90,) T laminate; a (90,/0,) T laminate; and a (0/90) 4S laminate. The first two lami- 
nates exhibit extremes in asymmetry, while the third laminate is a simple symmetric laminate 
with the same thickness and the same number of layers with fibers in both the axial and 
circumferential directions as the two unsymmetric laminates. 
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GEOMETRY AND COORDINATE SYSTEM 


The cylinder is assumed to be oriented in a global rectangular coordinate system with the 
X axis coincident with the centerline of the cylinder, as illustrated in fig. 1. The global Z axis 
is up, and the global Y axis is to the right. The origin of the coordinate system can be at the 
midspan of the cylinder, or at one end, depending on which location is convenient for the 
particular analysis. Here it shall be at midspan. The cylindrical coordinates used in the 
analysis consist of the axial coordinate, x, which is coincident with the X axis, 6 , which is 
measured positive from the +Z axis toward the +Y axis, and r , which is measured outward 
from the X axis. The cylinder has mean radius R, measured to midwall, and thickness H. The 
coordinate z is measured outward from the mean radius. The displacement in the axial di- 
rection is denoted u(x, 0,r), that in the circumferential direction (positive in the direction of 
+ 6 ) as v(x, 0,r), and that in the radial direction (positive outward) as w(x, 0, r) . The temper- 
ature change considered is assumed to be spatially uniform and is denoted as AT, AT being 
positive for temperature increases. Though interest here will be with applied axial end loads, 
the equations will be derived 'or the case of applied end torsional and applied radial loads. 


DERIVATION OF EQUILIBRIUM AND STABILITY EQUATIONS 

The total potential energy of a cylinder deformed by known applied loads, and subjected 
to preloading deformations, is written as 


7t = 4-|* [ [ ((<T X — "*■ { a 0 ~ °e) e 9 + ( T xtf — T x#)Vxtf) r d0drdx + TT'Qac (1) 

^ **r Jfl J x 

where 7 r lMd is the potential of the applied load. The stress components superscripted with a 
"P" denote preloading effects. They could be due to imperfections in the cylinder geometry, 
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for example, or due to thermally-induced deformations, as another example. They will be 
defined shortly. Using Donnell's assumptions for the kinematics of deformation, 

u(x, 0, r) = u°(x, 0) + z/?°(x, 0) 

v(x, 0, r) = v°(x, 0) + zpg(x, 0) (2) 

w(x, 0, r) = w°(x, 0) . 

In the above, as mentioned earlier, 


z = r - R 

and the superscript zero, as usual, signifies the displacements of the cylinder's reference 
surface (i.e., the surface at mean radius R), and the P's are the reference surface rotations 
given by 



0W° 

dx 

dw° 

R30 


(3) 


It should be noted that the displacements u°, v°, and w° are measured relative to the perfectly 
round cylinder, not the initial shape of the cylinder due to any preloading effects. The perti- 
nent strain-displacement relations in polar coordinates are 


_ __0U_ J_ / 0W \2 
dx + 2 \ dx ) 



_ dv du d\N 

y* 0 dx rdd + rdd dx 


Substituting eqs. 2 and 3 into eq. 4 and considering only thin shells so that the approximation 
r = R can be made with sufficient accuracy, the strains become 
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( 5 ) 


with 


e x “ e x Z *x 
0.0 
% = H + ZK * 

Vx0 = Vx0 ZK x0 * 


au° 

<3x 


±/ 

2 


0 

3v° 

, wl 

, J_ a 0 ' 

E e = 

' Rdd 

+ R 

+ ~2 Po 

0 

leu 

< 

o 

1 5u ° 

i n°n° 

- 

dx 

T Rdd 

+ PxPe 

o 

_ dfi 

0 

O 

K x = 

dx 

» k 9 

Rdd 


o 

*X0 


df 9 dp 


dx Rd9 


(6) 


The stresses are given by the relations 


o x = Qii(c x — £ x ) + Qi2( £ 0 — $ie(yx0 Vxff) 

a e = $12( £ x — £ x) + $22^0 — £ e) + $26^x0 — Vxfl) 
T xtf = $ie( £ x — £ x) + $26^0 — £ ff) + Qeefyxfl — Vxtf) » 


(7) 


where cj, c£, and are the strains due to preloading effects. Equation 7 can be rewritten as 


°x = $n £ x + $i2 £ 0 + $167x0 — a * 
a e = $12 £ x + $22 £ 0 + $267x0 — °d 
T x0 = $16 £ x + $26 £ 0 + $667x0 — T x0 


( 8 ) 


where 


<7 X = (Qh£ x + $12 £ 0 + $167x0) 
a 0 = ($12 £ x + $22 £ 0 + $267x0) 
= (0,8*: + $26 £ 0 + $667x0) 


(9) 


p 

T X0 


If the preloading effects are due to thermally-induced deformations, for example, then 
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( 10 ) 


£ x ~ «xAT 
cj = otflAT 

yls = a x0 AJ 

where a„ a B , and a* are the coefficients of thermal expansion in the cylindrical coordinate 
system. If this is the case, then the stress-strain relation of eq. 7 can be written as 

<*x = Qii( c x — <*xAT) + Qi2( c ^ — ®^AT) + Qi6(yxi> — a x^AT) 

a 9 ~ Ql2( £ x — a xAT) + Q22( £ 0 — a pAT) + Q26 (Vx 0 ~ Ot^AT) (1 1) 

T xtf = Qie( c x a xAT) + Q2$( £ 0 — a^AT) + QgeO/x# — a x6 AT) , 


or 

a x — Ql1 £ x + Ql2 e 0 + Ql67x0 — a x 

^0 = Ql2 c x + Q22 c 0 + ^26 yxff~ a 0 (12) 

T x0 = Ql6 £ x 4" ^26^0 "F ^66^x0 T x0 • 

where 

°x — (Qll“x + Ql2 a 0 + Ql6 a x0)AT 

a 9 = (Ql2«X + ^22*0 + Q26 a x0)AT (13) 

T x0 = (Ql6 a x + Q26 a 0 + Q66 a x0)AT . 

In this situation the superscript "T* denotes the fact that the preloading effects are thermally- 
induced. The stresses a \ , a ] , t* would then have the physical interpretation of being the 
stresses at a point if the composite is fully constrained from any deformation. 

For the case of known axial loads applied at either end of the cylinder 


‘load 


= J N x (0)u°( - ± , ey<id - J n£(0)u°( -t. , e^Rdd . 


(14) 
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The quantity N;(0) is the load at x = - y and N x + (0) is the load at x = + y. The integrals are 
taken around ends of the cylinder and the dimension of the N/s is force per unit 
circumferential length. For the case of a known torsional load applied at either end of the 
cylinder, 


*load - f NiWv°( - -t- , «)Rd9 - J N* (0)v°( (15) 
9 9 

and for a known outward radial load 


*ioad = - JJ" q( x - ^)w°(x, 0)Rd0dx. 


( 16 ) 


In general, the end loads can be functions of 0 and the radial load can be a function of 
xand0 The N*'s are tangential load per unit circumferential length, and q(x, 0) is a load per 
unit area. Other loadings on the ends of the cylinder can be included, e.g., an applied mo- 
ment, but they will not be considered here. 


Substituting the expressions for the strains, eq. 5, into the energy expressions for the 
cylinder, eq. 1, and including the three loading terms being considered, results in 


7r(u°, V°, w°) = — [ J J ((<7 X - <7 X )(£° + ZK°) + ((J ff - + ZKg) 

z "9 x 

+ ( T X* - ■&)(>& + ZK^))Rd0dzdx 
+ f N-(«)u°( - -t- , e)KM - J^N^(0)u°( i . «)Rd9 

+ J I , e)Rd9 - J 9 Ni(e)v°( -5- , «)Rd« 

— [ f q(x, 0)w°(x, 0)Rd0dx . 

-Vx 

In the above use has been made of the facts that 


(17) 
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r = z + R^R 
dr = dz 


( 18 ) 


and the integration on r has been replaced with integration on z . Integrating on z leads to 


rt(u°, v°, w °) - 2 J J (( N x “ + ( N e ~ nJ)c£ + (N^ - N^)y^ + (M x - M x )k' 

+ (M 0 - m£)k£ + (M^ - M^) Rd0dx 

+ JV(0)u°( - , 0)Rd0 - J N+(0)u °( ± , e)Rd0 


(19) 


+ 


l N ^( 0 ) v °( ~ Y ’ 9 ) Rdd ~ J N ^ (0)v °( \ ■ 0 ) Rd0 


0 v - / j 9 

Ur- 


- j q(x f 0)w°(x, 0)Rd0dx 


The stress resultants in eq. 19 are defined as 


£ 


r h/2 

N x = I ^x^Z ** + ^ 12^0 + ^16^x5 + ®11^x + ^12*0 + ®16 k xj 0 N 

J -H/2 

r h/2 

N^= o e 6z = A 12 £ x + A 22 £0 + A 26 y^ + B 12 k x + B 22 k£ + B 26 k^ — n£ 

J -H/2 

PH/2 

^xfl 3 1 A^g£ x + A 2 g£^ + Aggy^ + B 16 k x + B 2 gK^ + BggK^ 

J -H/2 

rn/2 

M x = I za x dz — B^£ x + B 12 £0 + B 16 y^ + D^k x + — M x 

J -H/2 
f H /2 

I za^dz = B 12 e x + B 22 &0 + B 26 y^ + D 12 k x + D 22 k^ + D 26 k^ — 

J -H/2 
f H/2 

I ZTxfidz = B 16 £ X + B 2 g£0 4* B 66 y^ + D^k x + D 2 @k^ + , 

j u /o 


( 20 ) 


'—H/2 


H being the cylinder wall thickness. In the above 
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( 21 ) 


r— r— 

N x = I 2 °xdz = I 2 (Q 1l£x P + Q-|2 £ 0 + Ql6Vxtf)^ z 

J_ Jl J_ _H_ 

2 2 

f— f— _ _ 

Ng=J .JdZ = j (Ql2 c x + $22% + Q26yxfl)^ Z 

~~2 ~~2 

N^=J 2 h T x0 dz = J 2 r (Ql6 e x + $26^0 + Q66Vxe) dz 

~~T ~~2 

r— r— 

M x = f 2 cr p zdz = [ 2 (Quel + Q 12 eJ + Qiev!*) zdz 

J Jj_ J H 

” 2 2 

r— r— 

=J 2 a P $ Z dz = J 2 h (Q 12 Cx + 022^ + Q 26 y^) zdz 


M l 


H 


2 

H 


'x0= [ 2 T x0 zdz = [ (Ql6 c x + Q26 £ 0 + Q66yxtf) zdz • 
J H J H 

““a “2 


These expressions are the so-called equivalent preloading stress resultants. Again, if the 
preloading effects are therma ly induced, 


r— r— 

Nx = N x=[ 2 °1<*Z = I 2 (Qiia x + Qi 2 «« + Qi6ajATdz 

J H J H 

2 2 

f— f— - - 

N£ = Njs 2 crjdz = 2 (Qi 2 <X x + Q 22 a d + Q 26 a xtf)^"*" dz 

J H J H 

2 2 

r— r— 

= N^s f 2 r^dz = I 2 (Q 16 a x + Q 26 <x 0 + Q^ajAIdz 
J H H 

’"2 2 

r— r— 

= M x = I 2 <r x zdz = I 2 (Q„a x + Qi6 a e + Qi 2 a^ Tzdz 
J _H_ J _H_ 

2 2 

r— p— 

= 2 a d zdz= 2 (Q 12 a x + Q 22 a tf + Q 2 6 a xfl) ATzdz 

J H J H 

2 2 

r— r— 

= M^= I 2 t^zdz = [ 2 (Q 16 a x + Qj6 a tf + Q66 a xo) ATzdz ■ 

J H J H 


M 


M 


( 22 ) 
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In this case these expressions are the so-called equivalent thermal stress resultants. 


The notation 


n = 7r(u°, v°, w°) (23) 

is being used to emphasize the fact that the total potential energy is a function of the dis- 
placements (here the cylinder reference surface displacements). The governing conditions 
will be derived by examining variations in the total potential energy. These variations in the 
total potential energy will be due to variations in these displacements. To this end, consider 
the increment, or variation, of the total potential energy due to increments, or variations, in the 
displacements. Specifically, 


u° + cu° ; v° + ev° ; w° + £w° , (24) 

where £ is a small parameter and the quantities u? , v? , and w? satisfy all the kinematic re- 
quirements of the problem. Using the notation 

71 + An = 7r(u° + £U°, V° + £V°, W° + £W°) , (25) 

the incremented total potential can be expanded using eq. 19 as follows: 
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n + An = y f J {(N x + AN X - N x )(e x + Ae°) + (N^ + AN# — N^)(c^ + Ae e ) 

+ (N^ + AN X<? - N^)(y^ + Ay^) + (M x + AM X - M x )(k x + Ak x ) 

+ (M 0 + AM* - mJ)(k£ + A 4) + (M* + AM* - M*)(k* + A»&)}Rd0dx 

+ f N-(S)ju°( . «) + tu?( - j- , e)}Rd9 

- i- . 6) + «i?( -t- , S)}Rd9 

+ J Ni(fl)|v°( - y •«) + «?(- -j- ■ e )} Rd « 

- J N i ( « ) {>'°( y ■ e ) + «?' ( T ■ ®)} Rd « 


- J J q(x, 0){w°(x, 0) + £w°(x, 0)}Rd0dx 

* X 0 


( 26 ) 


It should be noted that the equivalent preloading stress resultants do not have increments 
because they depend only on material properties and initial displacements (or perhaps the 
temperature), not the displacements due to the applied forces. Subtracting eq. 19 from eq. 
26 leads to an expression for the increment in the total potential energy, namely 


An — g J* J* {K — N <)Ae x 4- AN x c x + AN x Ac x 

+ (N* — N*)Ae* + AN*£* + AN*Ae* + (N^ N^^Ay^ "h AN^Ay^ 

+ (M x - M x )Ak x + AM x k x + AM x Ak x + (M* - m£)Ak£ + AM*k* + AM*Ak£ 
•+■ (Mjj* — M^^Ak^ + AM^k^ + AMj^Ak^I Rd0dx 


|J^n;(0)u?( - j- , o)Rd e - JV(0)u°( Y . 0)Rd0 J 

|[n^K( - j- , d)Rd0 - J(0)v?( y . *)Rd0 j 

- £ |j J q(x, 0)w°(x, 0)Rd0dx| . 


+ £ 


+ € 


(27) 
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The various increments in eq. 27 are given by substituting the increments in the dis- 
placements, eq. 24, into the basic definitions of strain and curvature, eq. 6. That substitution 
leads to, for e;, 


4 + = 




Dropping the superscript for convenience, 


+Al * = 'ir +£ ir + t^ +c/, a+-2- £ X 


Then, using eq. 6, 


^x = £ ( 5x 


5u 1 


+ PxP: 




(28) 


(29) 


(30) 


This can be written as 


— ££ x1 



(31) 


where and i s are defined to be 


3u 1 


£x i “ ax 






(32) 


(33) 


In a similar fashion, 


. 5(v + £V<) W + £W.< 1 , , 

c e + A c 0 = — + - + — (p e + 


(34) 


dv £vi 
R50 C R50 


W W-« 1 o 1 

f "r~ + £ -R - + Y Po + tPoPoi + ~2 



(35) 
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This results in 



C*5 

CO 


Cn 

rs» 


II 



X) Cl> 
Ql> < 
Q> -> 


+ 





+ 




to 


C*5 


M 




where 


Kv = 


Oft* 


* dx ’ 


(45) 


Kq ' RdQ ’ 


(46) 


d K 


i \ dx Rdd 


K X0 - 


)■ 


(47) 


In the above use has been made of the definitions 


A.- 


dw., 

3x 


! @$i ~ 


dw., 

R50 


(48) 


In terms of the strain increments, the increments in the stress resultants are: 

AN X = A^Acx 4- A 12 Ac 0 4 A^Ay^ 4 Ak x 4 B 12 Ak s 4 B^Ak^ 
AN^ = A 12 Ac x 4 A 22 Ac e 4 A 26 Ay x g 4 B 12 Ak x + B 22 Ak^ 4 B 26 Ak x5 
AN*# = ^ie^x + A 26 Ac fl + A 66 A y^g + B 16 Ak x 4 B 26 Ak 0 4 BggAjc,^ 
AM X = B^Ac x 4 B 12 Ae ? 4 B 16 Ay x ^ 4 D^Ak x 4 D 12 Ak^ 4 D-igA/c,^ 
AM^ = B 12 Ae x 4 B 22 Ae^ 4 B 26 Ay x5 4 D 12 A/c x 4 D 22 Ak 0 4 D 26 Ak x ^ 
AMxff = ®16^®X B 2 gA£^ 4* ®66^Vx0 ^16 ^ k x 4" D 2 g A/C^ 4 DggAK^ 


It is convenient to expand the increments in the stress resultants, and redefine those incre- 
ments in terms of powers of £. From eq. 49, incorporating the definitions for the strain incre- 
ments, eqs. 28-47, the stress resultant increments are 


AN X — A^ec^ 4 c fixj) 4 A^ce^ 4 £ £# 2 ) 4 A^fiy^ 4 £ y^) 
4 B^ck^ 4 B 12 ek^ i 4 B 16 ey x ^ i . 


(50a) 


These terms can be redefined to give 

AN X = eN Xi 4 £ , 


(50b) 
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where 


N Xi = A^e^ + A 12 c ffl + A 16 y xffi + 4- B 12 k Si 4- B 16 k x5i 

and 

N x* - A 11 c Xj + A 12 c 0j + A 16yxff 2 • 

Expanding for future reference, 

( du< \ ( w< \ 

Nx, = Al1 ( ~eT + ) + Al2 \ Rdd + T + ) 

+ A,e ( IT + 1m + fA < + w *-) 

a/L sh, ( «/>*, 

+ Bl1 ax + Bl2 Raa + Blf \ ax + Raa 

and 

Nx* = ( * 2 * A ii/*x, + * 2 * A i2^ + A-ieAx,/^,) • 

Using this procedure for the remaining stress resultants, 

AN* = A 12 (e£ Xi 4- £%) 4- A 22 (ce^ 4- e\) + A 26 (e yxffi + ty J 
4- B^ck^ 4- B 22 ek^ i + BjeEK^ . 


Redefining, 

AN* = eN^ 4- £ 2 N* 2 , 


with 


(50c) 


(50d) 


(50e) 


(500 


(51a) 


(51b) 


N^ = A 12 £ x, + A 22 £ 0, + A 26Vxd, + B 12 k x 1 + B 22 K 0, + ^26^x5, 


(51C) 



and 


+ A 22 £^ 4- A 26 y^ 


Expanding for future reference, 


and 


N«, = A, 2 | 

( du-, 
\ dx 

+ PxPx, 

+ Ajg( 

f dv, 

du* 

1 - ' 

^ dx 

Rdd 

+ B 12 

*K 

< 

I Q _ 

dx 

H- D 22 | 


^) +A “(w +J r + w«.) 


We y ( dfa, dp 


*26^ 


4 B?fil V h 


Rdd • ~ Z6 \ dx r Rdd 


and 


Likewise, 


or 


where 


( 2 ^12^X1 + 2 ^22^1 + ^2eP^0^j • 


AN^ A^ec^ + £ c Xj ) + A 26 (cc^ i 4- £ e^) 4- A 66 (£y x ^ i + ^YxsJ 
+ BigCK^ + B 2 g£K^ i 4* BggEK^ , 


4- £ N xfljt 


~ AigCx, + A 2 g£g i 4- Ag 8 y xtfi 4- B^k^ 4- B 26 K tfi 4- B 66 k x5i 


~ AigC^ + A 26 £ 5j + Agg y^ . 


(51 d) 


(51e) 


(5 If) 


(52a) 


(52b) 


(52c) 


(52d) 
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Expanding, 


N X0, - A 16 


+ A 


6 f! 


+ B 


If) 


dx 

5v 1 

dx 

dx 


\ ( dv* w., ' 

+ PJx,) + A 2 6^ -£qq + — + Me, 


du a 

TdQ + ^ 9 1 + 


(52e) 


+ B. 


26 


We, 

R6d 


+ B 


We, *K 


661 


dx 


R dd 


and 


N 





(52d) 


The increments in the moments can be similiarly defined, namely. 


AM X = + c 2 £xj) + B 12 (a #i + ti B ) + B 16 (£y^ t + ty^) 

+ Dh£K Xi + D 12 £K &1 + D 16 £K x( , i , 


(53a) 


or 


AM X = eM^ + e 2 !^ 


with 


= B 11 e Xl + B 12 £^ + B 16 y x ^ i + + D 12 »c tfl 4- D 16 k: x ^ i 


and 


+ B^ 2 £^ ? + BieVxtfj • 


where 


(53b) 


(53c) 


(53d) 


19 



(53e) 


Mx, — Bn 


Sua 


+ B 


161 


dx 

5v 1 


+ D 


11 


dx 

dx 


+ ^) +8 ’{w + ' 

5ui \ 

+ W + W «. + Wx,j 


W, 


+ PoPo, 


dp 9 i 

+ D,J 'Ra5‘ +D ’ 6 


«/>., <%< 


dx 


d K \ 

Rdd ) 


and 

Mx, = ( ~2 Bii^x, + ' 2 ’ ®i2^, + • 

In a similar manner, 

AMp = B 12 (ec Xl + e £ xj) + B 2 2( ££ 0 1 + £ + B^cy*^ + £ y^) 

+ D^ck^ + D 22 £K^ f + D 2 g£K x ^^ , 


or 


AM tf = + £ 2 M^ , 


with 


— Bi 2 £ Xi + B22 e 0 1 + B 26 y xfli + D 12 k Xi + D 22 k^ + D 26 k x ^ i 


and 


— Bi 2 £xj + B 22 £ (9j + B 26 y xtfi , 


where 


(530 


(54a) 


(54b) 


(54c) 


(54d) 
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and 


Finally 


where 


du< 

= B12I 


) + B22 ( Rdl 


Rdd R 


+ ^~ + hh, 


+ B26 ( dx + Rdd + W 0 ' + ^ x \ 


+ D 




We, ( We, W* 


12 


F D 2 2 TTTX" + D 


ax 1 22 Rde 


261 


ax 1 Raa 


(54e) 


M, 


= ( \ B„/g, + \ B J2 4, + 


(540 


AM^ = B^ee^ + c 2 ej + B 26 (cc 01 + £%J + B^Ey^ + e^J 

+ D^CK^ + 026^^, + 


(55a) 


AM^ = eM^ + c M^ z 


(55b) 


with 


M xtfi = B^^ + B 26 £ Si + B 66 y x ^ i + Di6 K x, + D26*#, + D 66 k x5i (55c) 


and 


M xfl 2 = B 16S + B 26 £ <? 2 + B 66Vxtf 2 ■ 


(55d) 


Expanding for future reference 
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diii 


Pxfi) 


^ + B 2e(‘ 




Wi 


8x • r^j ■ -wy Rgd 

_ / 5v 1 du-. \ 

+ Be6 V "dx” + W + + ^ x lJ 


(55e) 


+ D 


16 


fAj, 

dx 


+ D 


26 


Rdd 


+ D fi 


*K W* 


dx 


Rdd 


and 


M 


X0 2 



(550 


With the various increments defined and expanded, the definitions can be substituted into eq. 
27. This results in 


= T J J {(Nx - <)(««„, + e\) + 

( £ ^x, + c N Xj )(eC Xi + £ Cxj) + (N0 — N,)^^ + £ 2 £^) 

+ (eN 9i + £ 2 N 0j )£ 9 + (fiN fli + £ 2 N <? J(££ (?i + £%J 
+ (N^ - N^,)(£y x , 1 + e 2 Vx« 2 ) + ( £ N x , 1 4- £ 2 N x , j )y x , 

+ (eN x , i + e 2 N x , 2 )(£y x , i + t 2 y^ + (M x - 
c Nl^Kx 4- (cM^ 4- £ MxJeKx, 

+ (M, - M %) ck 9i + (eM^ 4- c 2 M, 2 )k, + (eM^ + £ 2 M, 2 )£K, f 
+ (Mx, - M^)£K X , 1 4 (eM x , i 4 £ 2 Mx, 2 )>Cx, + (cM^ 4- ^M^JeK^jRd^dx 

+ £ |j^N; (0 )u?( - ^ , O^Rdd - |n+(0)u ?( ± , d^Rdd J 

+ £ jjj N*(0)v?( - y , d^Rdd - |n+ (0)v?( \ , 0)RdQ 1 


e|J J q(x, 0) w°(x, 0)Rd0dx j . 


) 


(56) 


Expanding and regrouping in powers of £ leads to 


22 


Art = y f J [e{(N x - + N^Ex + (N^ - N^e^ 4- N^e# 

+ (N* - iy x ^ i + + (M x - Mx)^ + M^Kx 

+ (M tf - Mj)*^ + M e Kg + (M x0 - M^)k x9i + M^k^} 

+ e 2 {(N x - N^Ex, + N X2 e x + N Xi e Xi + (N tf - nJ)e* 2 + N 9 t g + N \ 9 fa 
+ (Nx, - N^ ly^ + N^y^ + N^y^, + M xj*x + M x, K x, 

+ M 0j K e + + M^Kxtf + 

+ c 3 {N Xi£xj + N^Ex, + N 8 t 0t + + N + N^y^ 

+ M Xj >Cx 1 + + £ {Nx^^ + N0 2 c 0j 

+ Nx^jjRdfldx 

+ e |j n;(9)u?( - i. , #)r<» - JViJ(e)u?( \ . e)Rd#| 

+ « jj < Ni(9)v?( -±,e)M9- J # Ni(9)v?( ±. . »)Rd«| ■ 

- e{ f f q(x, 0)w?(x, 0)Rd0dxl . 


(57) 


The increment in the total potential energy can be written as 


2 3 4 

Art = E7T, + E rt 2 + E rt 3 + E rt 4 


(58) 


The quantities n v n 3 , and rt 4 are defined to be the first, second, third, and fourth variation, 
respectively. The equilibrium conditions for the cylinder are obtained from the condition 

rt^u,, v 1t w.,) = 0 , (59) 


where the notation indicates jr, is to be made stationary with respect to the displacements 
Ul , v„ w v These displacements are the variations in the equilibrium displacements. The sec- 
ond variation is used to examine stability of the equilibrium displacements. According to the 
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Trefftz stability criterion, transition from a stable equilibrium configuration to an unstable one 
is characterized by 


^(Uv^.wO-O. (60) 

This states the second variation of the total potential energy should be stationary with respect 
to variations in u,, v„ w,. The above two conditions will now be examined. 

First Variation 

The first variation can be isolated from An and is 


*1 2 J J ^x, £ x + (N# — 

* x ^ 0 

+ (N* - + N^y* + (M x - 

+ (M, - m£)k 9i 4- M 0 k 9 + (M a + M^/c^Rd^dx 



(61) 


A more useful form of the first variation can be obtained by substituting for 
N*, . N #1 , N rfi , , M #1 , and M rf , from eqs. 50e, 51e, 52e, 53e, 54e, and 55e. If this is done 

and the various terms in this expanded form of n y regrouped, the result is 
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( 62 ) 


*i=i| {N^ + N,e, + + M x k Xi 

J x J 0 

+ M qKq^ — M^K^jRdfldx 

+ | |n-(8)u?( — |- , «)Rd» - J # Nx(«)u?( y • »)Rdff| 
+ jjki(»)v?( — . e)Rd» - Jn J(0)v?( \ . e)Rd«| 
q(x, 0)w?(x, 0)Rd0dx > . 

x-'fl J 



Note the quantities N* have disappeared, as has the factor of 1/2. 


If the strain and curvature increments t ^ ,... k are written in terms of the displacement 
increments, using eqs. 32, 37, 42, and 45-47, the first variation takes the form 




6v 1 

+— + 


+ N **( IT + W + W*' + + 

dfai dP* 


* Rdd ' r 


■) 


M, 


dx 


+ M* 


Rdd 


+ M 


X0\ 


dx Rd6 


>Rd0dx 


(63) 


+ JJn'u, ( — t- , e) - nJu,( \ , #)}Rd« 

+ J |n>i( - y . «) - N >i( y • 9 )} Rd(i - J |qMw,(x, SjRdfld) 


This is one of the fundamental forms of the first variation for determining the response of a 
cylinder. This form can be used directly in approximate schemes such as the Rayieigh-Ritz 
method. However, here we are interested in the governing equilibrium equations and asso- 
ciated boundary conditions. The steps to derive these follow: 
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To determine the equilibrium equations, and the associated boundary conditions, differ- 
entiation of u 1 ,v 1 , and w, with respect to the spatial variables x and 6 must eliminated. This 
is done using integration by parts on the various terms in eq. 63. This shall be done in the 
following term by term: 

first term 



du 1 

N x — — Rd0dx 
dx 


- f (N„u.)| 

J 9 

dH 


X “ + T 


Rd0 


-LI 


-r — u,Rd0dx 
dx 1 


(64) 


second term 


1 A Rdf,dx - £ J s N x 3 d “ 

=I( N ‘17 w ')Cl Rde 

“JJ # ir( N »- fr) w ’ Rd ' Wx - 

In the above use has been made of eqs. 3 and 48. 
third term 



Rd0dx 


-L 

-II 


(N^) 

5N, 


| 0=+fT 


dx 


It 


w v,Rd8dx. 


( 66 ) 
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Because the cylinder is complete, the response is 3 continuous function at the spatial variable 


0. Thus 


Nfl(x, + n) = N^(x, - tt) 

v^x, + n) = v^x, - 7t) . 


(67) 


As a result the first integral on x in eq. 66 sums to zero and we are left with the third term as 



Rd0dx = — 



R00 


v n Rd0dx . 


( 68 ) 


fifth term 


f f 

** 0 


Rd0dx 


■U" 

-Ji 


0w d w i 
6 R 69 R00 


Rd0dx 


R00 




Rd0dx 


(69) 


where use has been made of eqs. 3 and 48, and the fact that the first terms normally on the 
right when using integration by parts is zero because the cylinder is complete. 


sixth term 




f f 1!1 *l 

"J x J e 0x 


(70) 
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seventh term 


eighth term 


ninth term 


l l Nx » Tde Rd0dx = " J x // Tw - u ’ Rd0dx 


(71) 


l l N **M*i Rd0dx - 1 l N x* ax ~Rdd Rd0dx 

=_ !J 0 w^'lr ) w i Rd0dx 


(72) 


f f N x ,/?,/? Xi Rd0dx= f f N* 
J x J e J x J e 


3w ^w 1 

rbq ax 


Rd0dx 


-JliK"-™)- 


Rd0 


Rdadx 


(73) 


tenth term 



Rd0dx = 




aM, 

dx 



Rd^dx 


(74) 
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These terms can be expanded further. Specifically, using eg. 48, the second term on the right 
side of eq. 74 becomes 




5M X dw., 
dx dx 


Rd0dx 


(75) 


Using integration by parts once more yields 


-[/■£ 

V* dx 


dx ^ 


/LRdfldx 


-II 


W ’Jl” + i RdS 

2 

3 2 M, 


w.,Rd0dx 


(76) 


The tenth term can thus be written as 


J x J e 


dx 


Rd0dx 


-I,("'Tr)Ci 


Rd0 


3M X \ix=+t 
- r-^- w< 2 Rd6 

V'x-i. 


-u 


a 2 M, 

ax 2 


w^dddx 


(77) 


eleventh term 



Rd0dx = — 



d M, 


W*. RdMx ' 


(78) 


where continuity of the cylinder has been used. The expression on the right can be further 
integrated by using the definition of /? #1 from eq. 48, i.e.. 
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n dM 0 n f r dMg 5w, 

. J wV dte= |J,Wi» Mdx 


d 2 U 


9 


— 7 — — w.,Rd 0 dx 
j R 2 dd 2 


( 79 ) 


Thus the eleventh term becomes 


twelfth term 


a 


Wt, 


M »w RdWx 


-JJ 

Ja 


d 2 M 


9 


, , w«Rd0dx 
x-V R 2 dd 2 


(80) 


j x J^-^Rd0dx = ^(M^ i )|^| Rd0 


-ffi^ 

J X J 9 dx 


Rd0dx 

dx 


(81) 


The second term on the right can be further integrated as 


n dM *0 „ r r <9M xfl a Wl 

. ^ Rdedx = lJ # -5f W Rdedx 


-u 


R00dx 


(82) 


w.,Rd0dx 


where again completeness of the cylinder has been used. Hence the twelfth term becomes 



-f(M^rr Rd0 

J 9 x =-y 


n d 2 M xff 

. J Mwr w ' Rdddx - 


(83) 
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The first term on the right in this resulting equation is special. It can be integrated by parts 
one more time. Considering the integral at x = + L/2, the term becomes 


J M x# ( + ± ■ «K( + T • s ) Rde “ - J 1 + T’ 9 ) 


<9w 


•(♦-M 


R89 


Rdd 


(84) 


Integrating by parts 


J M*( + \ .«)*,(+ Y • e)R d » -- m*( + j - , + *) 

e 

+ Mx*( + y . - *) - 


5Wl ( + "2*’ + *) 


Rdd 




R89 


C 3M*( + A- . #) 

+ J tf Rdd 


+ -^-,0)Rd0 


(85) 


Because of completeness of the cylinder, the first two terms on the right sum to zero. The 
same procedure can be used at x = — L/2. The result is 



Thus the final form for the twelfth term becomes 


(86) 



17 d M x<? \ x-+y 
“U Ra0 7 x=-i 

- f f w 1 Rdfldx . 

JJa R503x 1 


(87) 
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thirteenth term 


JJ 9 M "'^- Rd9dx = -JJ # W^ Rdedx ’ < 88 > 

where completeness of the cylinder has been used once more. The term on the right can be 
further integrated by the fact that 


n dM 

j ™ 


dM„ 0 dw« 


^ Rd8dx - 1 1 Rd8dx 
■I( “ 


-)rt 1 

/ X— T 

, Rdddx w i Rd0dx 


Rdd 

d 2 M 


(89) 


The thirteenth term thus becomes 



Rd0dx = 



dM„a 
~Rdd Wl 


Rd0dx 



w.,Rd0dx 


Rd0 


(90) 


Using eqs. 64, 65, 68, 69, 70, 71, 72, 73, 77, 86, 87, and 90 in the expression for the first vari 
ation, eq. 63, and combining boundary terms (i.e., integrals on 6), results in 
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For the first variation to be zero, each individual term in each of the integrals must be zero, 
the Euler equilibrium equations coming from the two-dimensional integral, and the boundary 
conditions coming from the one dimensional integral with respect to 0. Thus, the three gov- 
erning equilibrium equations are 


0N X 0N X(? 
dx + Rdd 


(92a) 


dN X 0 

dx 


SN, 


Rdd 


(92b) 
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(92c) 


d i) 2 M x 

dx 2 


+ 2 


d 2 M 


xO 


d 2 M f 


+ 


Rdddx r 2 50 2 


Rdd 

Ni 


a 

R50 


4-q = 0 



N 


x0 


aw \ 
Rd6 ) 


Using the first two equilibrium equations in the third one, the three equations can be written 
as 


dN x aiSL„ 

~ + w =0 ( S3a ) 


dN M <?n„ 

— + W = ° (93b) 


£m, 

ox 2 


4-2 


5 2 M 


x<9 


+ N, 


Raaax 

a 2 w 

0 R 2 aa 2 


_a^M 

R 2 aa 


r+N,-^f 

2 ax 2 


+ 2N 


xB 


N* 

if + q=° 


a 2 w 

Raaax 


(93c) 


The variationally consistent boundary conditions at the ends of the cylinder are: 
at x = - L/2 

i) N x = N x or u must be specified , 

ii) N„0 = or v must be specified , 

*'') 3x + "a7 + N ** ~Rdd + ^ " Ra^ = ® or w must be specified , ( 94a ) 

iv) M x = 0 or ~L must be specified . 
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at x - + L/2 


i) n x = N x or u must be specified , 

ii) or v must be specified , 

_5w_ N _5w_ 2 d^xg = o or w mu st be specified , 

ax x dx * Rdd Rdd 

iv) M v = 0 or ~~ must be specified . 

' x dx 


(94b) 


Second Variation 


With the equilibrium conditions established, attention now turns to the stability of these 
conditions. As stated by the Trefftz criterion, stability information can be obtained by exam- 
ining the first variation of the second variation, n v The second variation can be isolated from 
A n of eq. 57 and is given by 


n 2 = \ J J [( N x - <) Cx2 + fV x + + (N, - NJK 

+ IV* + I V*, + ( N xg ~ *&) V*, + N xg 2 V^ + (95) 

-I- + M 6 k 9 + M 0 k 0 , + 

+ M^K^jRdfldx . 

A more useful form of n t can be obtained by substituting for 

N,,. N tfj , N* 2 , M, f M, 2 ,and from eqs. 50f through 55f. Doing this, and regrouping 

terms, results in 


J x J 0 

+ J (N^x, + N + Mx,*x, + M 
+ Mxg «xgj] Rd0dx 


(96) 
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Again note the disappearance of NJ ,NJ , and N£, and the algebra with the factor of 1/2. Since 
the variational process involves kinematic variables, here u, ,v, , and w„ the increments in the 
strains, curvatures, and stress resultants should be written in terms of these kinematic vari- 
ables before the variational process begins. To that end, using the expressions for 

N*i >—* eqs. 50e through 55e, k ^ from eqs. 32, 37, 42, and 45-47, and 

. and from eqs. 33, 38, and 43, 
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Using the definition of N v N Sl M -1t and combining terms leads to an important form of eq. 


98, namely. 


<5tt 2 = J J + N ePgSfie, + N + Pe&PJ + n *i('1T' + 

+ N(? i( 1 m~ + _ V L + p ° S(i9 ') + Nxtf, ( ~dT + W + P* 5 ^ 9 ' + Pe5fi *' 


+ M 


R50 

M0*, 


(99) 


* dx 


s ''Rar +M "| _ Ar + R3» )} 9 


Since stability is studied by the condition <57r 2 = 0, the above integral is equated zero. Then 
this form of <57r 2 = 0 can be used for approximate solutions to the stability conditions. A 
Rayleigh-Ritz formulation can start from this form of Sn t = 0. 

To determine the differential equations, and associated boundary conditions from which 
to study stability, integration by parts is used to eliminate differentiation of the variables of 
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Ut,v 1( and w, with respect to the spatial variables x and d This is done, as before, on a term 
by term basis as follows: 

first term 


l l N A M*, RMdx - £ J # N » ^ ax' Rd<Wx 

J x «w,) | ~ r RM 

-J x iir( N *'ir)' 5w ’ Rd(,dx 

second term 

Rdfldx (102) 

l J, N * V % > R<Wdx = - { J # w ( N « -555- }*».> ROSdx . ( 1 03 ) 

where continuity of the cylinder has been used. 
third term 

ll N *Wh’ RMdX = U s N «^l^ Rdt>d * (104) 


( 100 ) 


( 101 ) 
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_[ j > N*V/l,,Rd«dx - - |J s W (Nx, )«w,Rd«dx 


(105) 


fourth term 


| J^/J^Rdddx-JJ^N 


aw., 55w 1 
5x“ 


Rd0dx 


(106) 


( J # N »>* ^«. RdSdx - J 9 ( N *» w 5w ’) I Z- 1 


Rd0 




(107) 


fifth term 


U N «.-ir Rd< ’ dx -J» (N ^ u ' ) 




<5u,Rd0dx 


(108) 


sixth term 


JJfN^WdWx-^Jfn 


aw ^w, 

*1 0x dx 


Rd0dx 


(109) 


41 


seventh term 


ninth term 


JJ^xiW^RdWx = J (n„ ■§*-*,,) RdS 

2 

~JJ. ir(^lr} i,v ’ RdMx 


(110) 


n ddv^ r r dN ff 

. >wr Rd 9 dx= -JJ # w 5v ’ Rd9dx 


(111) 


f f \M0e. Rd^dx =[ [ N« 

J x J 0 J x J 0 


Rdd Rdd 


Rd0dx 


(112) 


I x J # N * AM>, Rd « dx - _[ J # (n„, )**,RdSdx 


(113) 


tenth term 



aav 1 

dx 


Rd0dx = 



Rd0 



<5v.,Rd0dx 


( 114 ) 
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eleventh term 


JK 

J x J e 


ddu^ 

Rdd 


Rdfldx 


__f f aN >*i 

_ J X J* Rdd 


<5u.,Rd0dx 


(115) 


twelfth term 


J | N xS ,W/!,,Rd0dx = JJ < N x# , 


dw ^ w i 
dx Rdd 


Rd0dx 


(116) 


( j # N^w/f.Rdedx = - J x J # w ( N >», 17 )‘ 5w ' Rdedx (1 1 7) 


thirteenth term 


J J N^/V^RdCdx = 


dw ^ w i 
R00 dx 


Rd0dx 


(118) 


f ( IVWdOdX - f («*, -gjg ^ 1 

J X J $ ' 


x=+-^ 


5F aw ’;i„. + 


Rd0 


-U^r( N -.^) iw ’ Rd9dx 


(119) 


fourteenth term 



£^x, 

dx 


Rdddx = 





SPxRdddx 


( 120 ) 
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Using the definition of /? v and interchanging the <5 and differential operators, the second term 
on the right hand side becomes 


-u 


<5Mx 

dx 


■d/LRdddx 


j dx 


dS w, 
dx dx 


Rdddx 


( 121 ) 


or 


-u 


d _ Mx 
dx 


■ d/LRdddx 


■/.( 


<?M X , 

dx 


<5w 1 


-i 



d 2 M x 

ir-dw^dddx 

dx 2 


Rdd 


(122) 


With this the fourteenth term becomes 



Mix, 

dx 


Rdddx = — 



Rdd 


+ 






d 2 M„ 

r— dw.,Rdddx 

dx 2 


(123) 


fifteenth term 


n dSp e ^ r <• 

. J M ».w Rd9dx= -JJ s w^ Rd( ' dx o»> 


Using the definition of /? #1 , and interchanging operators, the right hand side becomes 


-Uw^ Rdddx =JJ 


A 


(125) 
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or 


n SM 

,Td 


a 2 M„ 


wir ^ Rd( ’ dx " - U ^ {w ’ RdWx 


rw 


(126) 


Therefore, the fifteenth term becomes 


n ddBa f f ^ 

M 0 p Jr -Rdgdx = - g Wl Rd0dx 

„ j 6 ' R50 J X J 0 R 2 d0 2 


d 2 M fi 


(127) 


sixteenth term 


r r 

-U.sr-'W** 


(128) 


Both of these integrals can be expanded by using the definition /? #1 . The first integral becomes 


f x -+7r f/ 5< ’ w i \l X= * + '2' 

J (^W^) ± Rd 0 = -J^ M x<? 1 Rdd )\^__L 


Rd9 


(129) 


**■ 




’JL-l 


(130) 


where continuity of the cylinder has been used once again. The second integral on the right 
hand side of eq. 128 becomes 


U 5M 
3~‘ 




(131) 
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or 


-a 


— ^ Rd9dx —11 ' 5w ' Rded) ' 




(132) 


Hence the sixteenth term becomes 


l/>^—l(^-.)Ci Ms 

“11 M^r ■ 5w ’ Rdedx 


(133) 


seventeenth term 


Rewriting /? , 


IK 


dS(J Xl 

Rdd 


Rdddx 


-^-dfi^Rdddx 

* 0 


Rdd ^ 


(134) 


I K 


5<5 K f f aM x#, d<5w., 

Rd0dx = J J “5*2 ^Rdfldx, 


Rdd 


Rdd dx 


(135) 


the seventeenth term becomes 


Rd0 

(136) 

n ° Mx». 

. , Rae3T ,iw ' RdWx 

Using the results from eqs. 101, 103, 105, 107, 108, 110, 111, 113, 114, 115, 117, 119, 123, 127, 
133, and 136, the variation of n t becomes 
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Based on the above variation, the Euler equations are 


£Nxj_ dN*, 
dx + Rdd 


(138a) 


di\ 

dx + Rdd 


(138b) 


£X 

dx 2 


+ 2 Rdddx + r 2 ^0 2 


+ 


/ 5w 1 5w 1 5w 

ar( Nx ~ + "W +N,< ' ax + 

d ( 8 ' N ^ ^ w i , ki dw 

^f( N 9 "r^* + n ^ 8 * N<?i R50 


dx 

+ - Rd0 


'x0. 


+ N, 



(138c) 


Using the first two equilibrium equations from setting n y = 0, eqs. 93a and b, and the first two 
of the above in the third equation above, the three Euler equations from 3n t = 0 can be written 


as 
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(139a) 


5N„ 

21 

dx 


dN 


X0. 


Rdd 


= 0 


3N*, 8N e 

-H = n 

dx * Rdd u 


gX 

dx 2 


+ 2 


Rd0dx 


R 2 dd 2 


+ N X 

dx 2 


5 2 Wi 3 2 \n 1 


nil O w 1 d 2 w, 

2Nvfl -zrrr-z + N 0 - 

** Rdfldx + * R 2 a0 2 


, — 4- OKI w d 2 W 

5(1 5x 2 *’ R^x + Ntf i R 2 a0 2 


+ N^ 


N, 


= 0 


The variationally consistent boundary conditions are: 


(139b) 


(139c) 


at x = + L/2 


i) 

ii) 

iii) 


vi) 


= 0 or u., must be specified 
Nxfl, = 0 or v, must be specified 


(n, 


d w. 


+ N, 


d w 




+ N, 


d w 


* dx Rdd " ,,x i dx 

or w., must be specified 


+ N 


dw 
Rdd 


= 0 or 


dw, 

dx 


must be specified 


dx 


dM 


+ 2 - 


Xtf, 


Rdd 


= 0 


(140) 


Equations 139, or the alternative form eq. 138, are referred to as the buckling equations. They 
along with the boundary conditions of eq. 140 provide the conditions that must prevail when 
the cylinder passes from a stable equilibrium configuration to an unstable configuration. Sol- 
ution of eqs. 139 satisfying eqs. 140 leads to the value(s) of the applied load(s) that cause in- 
stability. Solution of these equations also gives the buckling shapes relative to the equilibrium 
configuration. 


48 


SIMPLIFICATION OF THE EQUILIBRIUM EQUATIONS DUE TO THE CONDITIONS 

OF AXISYMMETRY 

The studies here will focus on the case of axisymmetric end loads and axisymmetric re 
sponse. For this situation 


g ) 

dd 


= o 


and 


d( ) _ d( ) 
dx dx ’ 


(141) 


( ) being any response quantity. With this condition, the kinematic relations simplify con- 

siderably. Specifically, using eqs. 3 and 6, and reintroducing the superscript 'o' notation 
where applicable, 


*x = 


£ — 

1 


dw 

dx 


Pe=0 


dx 2 Px 


3 ° 

o 


dx 


c e - 


w 

R 


O dv 
y *° dx 


k° = 0; 




As a result of eq. 141, the equilibrium equations, eq. 93, simplify to 


dNx 

dx 


= 0 


dN x< > 

dx 


= 0 


d 2 M x 

dx 2 


+ 


N 


d 2 w 

X . 2 

dx 



q = 0. 


The accompanying boundary conditions are 


(142) 


(143a) 


(143b) 


(143c) 
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at x = - L/2 


i) N x = N x or u° must be specified , 

ii) or v° must be specified , 


iii) 


dM„ 


dw c 


dx dx 

dw° 


= 0 or w° must be specified , 


iv) M x = 0 or 


dx 


must be specified . 


(144a) 


at x = + L/2 


i) N x = N x or u° must be specified , 

ii) or v° must be specified , 


iii) 


dM 


x 1 1 dw o 


dx 


+ ISL 


dx 

dw' 


= 0 or w must be specified , 


iv) M x = 0 or ■ must be specified . 


(144b) 


Equations 143 and 144 will be the focus of the remainder of the study. In the next section the 
solution of these equations for the case of a known axial end load will be derived. 

SOLUTION OF EQUATIONS FOR THE CASE OF AN AXIAL END LOAD 


The first equilibrium equation of eq. 143 integrates to become 


N x = constant . 


(145) 


Since the axial load is known at the end of the cylinder, this constant is the applied end load. 
It will be referred to simply as N. The second equation of eq. 143 integrates to 


= another constant . 


(146) 


Since in the present work there will be no torsional loading on the end of the cylinder, the 
shear N* is zero there and thus this constant must be zero, i.e., 
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To solve the third equation of eq. 143 it is most convenient to express all quantities in that 
equation in terms of w° . This is accomplished by using the definitions of the stress resultants, 
eq. 20. Here the preloading condition is thermal. Thus eq. 22 is used in conjunction with eq. 
20. In using these equations it convenient to consider only cylinders that will not experience 
a global twist due to the application of axial end loads, i.e., consider the case of A„ and A a 
equal to zero. This is a valid assumption for most applications of cylinders. With this situ- 
ation, using the definitions of the stress resultants and the simplified kinematic relations, 


N x — A n e x + A 


w 
12 R 


- B 




dV -Nl-N 


dx' 


Nfl = A 12 C x + A 22 - B 12 


d 2 w° M r 
— N/> 


dx 


Nx# — A 66Vx0 B 16 2 

dx 


d 2 w° m t n 
-N x< , = 0 


(148) 


M x =B 11 c'< + + B 16 y^ — D n , M x . 


dx 


Solving the first equation for t\. 


and the third for yj,, 


A n 


N + N ' - A 


w 
12 R 


+ B 

+ " dx 2 )' 


(149a) 


o 


Bie dV 

A 66 dX 2 


N*L 

A 6’5 


(149b) 


Substituting these into the expressions for N # and M„ yields 


N* = 


'12 


'22 


W 


'11 


^7 (N + Nl) - Nj 


B 


12 


BjiA 12 

An 


d w 
dx 2 


(150a) 


M x = 


B 


BnA 


12 


11 ~12 
Si 


B 


+ 7^( N + N ^) “ Ml 4 



(150b) 


For the three laminates of interest, i.e., the (0,/90,) T , the (90,/0,) T> and the (0/90) 4S laminates, 


®22 ®ii > ®i2 Bis — B 26 — Bee = N x o = 0. 


(151) 


The expressions for N, and M, simplify to 


— ( A 22 — 


'12 \ w 


An 


+ ^7( n + n 3- n J 


A„B. 


1 2 D 11 

Si 


d 2 w° 

dx 2 


(152a) 


M x = - 


B 


A12B11 


w 


- D 


B 


11 


'11 


11 


+ ^L {n + n1) _ m t 


dV 


Si / dx" 


(152b) 


Substituting these into the third equilibrium equation, and considering the case of q = 0, leads 


to the equation governing w°. This equation is 




+ 


/ 2A 12 B 11 

V A n R 



d 2 w° 

dx 2 


Nj 

R 


A12 

RA n 


(N + Nj) . 


(153) 
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This equation is a linear differential equation with constant and known coefficients. The 
applied axial load N is known, as are the material properties, equivalent thermal loads, and 
geometry. This structurally nonlinear problem results in a mathematically linear problem. 


Having addressed the three equilibrium equations, attention focuses on the boundary 
conditions. Two of the four on each end have been specified. In particular, statements have 
been made regarding N x and N„. These satisfy eq. 144a and b, i and ii. To be considered are 
eq. 144a and b, iii and iv. Regarding these two remaining conditions, within the context of the 
admissible conditions, three physically plausible boundary conditions can be imposed on the 
ends of the cylinder. These are: 

1 - lubricated boundaries; 

2 - simply supported boundaries; and, 

3 - clamped boundaries. 

For lubricated boundaries the shear force and moment at the ends are zero. For simply 
supported boundaries the radial displacement and the moment are zero. For the clamped 
boundaries the radial displacement and the slope are zero. The terminology 'lubricated 
comes from the fact that the radial displacement has no resistance, as if a highly lubricated 
flat plate were pushing axially against the ends of the cylinder. This lubricated flat plate would 
not resist rotation of the ends either. 


For a lubricated boundary, the conditions are formally, from eq. 144, 


and 


dM x 

dx 


+ N 


dw° 

dx 


= 0 


M x = 0, 


(154) 


while for a simple support the conditions are 
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For a clamped support the conditions are 


w° = 0 


dw 

dx 


= 0 . 


(156) 


In terms of w°, for the three laminates of interest the lubricated support conditions are 


and 


111 


dw 

dx 


dV , ( A„B„ .A 
q _®11 \ d 2 w° / \ 


B 


- "T~ (N + Nx) + Mx = 0 . 
^ 1 1 


= 0 


(157) 


For a simply supported boundary, the conditions are: 


and 



(158) 


For a clamped boundary the conditions are as given in eq. 156. Here interest will focus on 
simply supported and clamped boundaries. 


PRELOADING RESPONSE DUE TO THERMAL EFFECTS 

Cylinders are generally fabricated on a mandrel and are consolidated at an elevated 
temperature. After consolidation the temperature is lowered to the ambient temperature and 


the finished cylinder is removed from the mandrel. Assuming that the fabrication and con- 
solidation are axisymmetric, the shape of the cylinder can be determined using the governing 
conditions from the previous section. Specifically, since there are no loads applied to a cyl- 
inder that has been removed from the mandrel, eq. 153 with N equated to zero governs the 
response. That is, 


D 




Bii \ d 4 w° 2A 12 B 1 i d 2 w° 

An J dx 4 AhR dx 2 



(159) 


With the cylinder simply in the ambient environment, the ends of the cylinder are traction free. 


Thus the two other boundary conditions, N, = 0 = N* being the first two, are 


dM x 

dx 


M x = 0. 


(160) 


In terms of displacements, the boundary conditions at x — ±L/2 are 



(161) 


The complete solution to eq. 159 consists of homogeneous and particular parts, i.e., 

w°(x) = w£ omo (x) + w part (x) . (162) 


By inspection, the particular solution is 
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(163) 


w part( x ) “ 


R(AjjN J- A 12 nJ) 


A 11 A 22 " A: 


12 


= W 


T 

part * 


where the notation wj. rt denotes the particular solution to the thermal problem. Note that it is 
not a function of x. 


The homogeneous solution is of the form 


^homoM — A e 


(164) 


Substituting this assumed form into eq. 159 results in the characteristic equation, namely 



AiiAga — A 
R 2 



= 0 . 


(165) 


The four roots to this equation are 


^1,2, 3, 4 - ±, 


A 12^11 \/ A 1 1 (D 1 i A 1 2 A 22^11 A 11 A 22^1l) 


^(DiiAii — B 2 .,) 


(166) 


These four roots are of the form 


^1,2, 3 , 4 = ± « ± ijff . 


(167) 


a and /? real and positive. It should be noted that the characteristic roots are strictly functions 
of the elastic properties and geometry. Thermal effects, in terms of material expansion coef- 
ficients or temperature, are not involved. At this point in the analysis a and /? are known. 

Because of eq. 167, the complete solution for the thermally-induced response is given by 
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w?(x) = A ie (0+l/?)x + A 2 e (a - Wx 4- A 3 e ( -“ +Wx 


+ A 4 e (_a_l/?)x + 


R (A-i i a 12 n x ) 


AhA 22 


-Ai 2 


(168) 


where A, are the four unknown constants of integration and the subscript T with w°(x) is used 
to denote the fact the solution is for the case of thermally-induced preloading. 

Since the boundary conditions, eq. 161, are the same at each end of the cylinder, the 
thermally-induced response is expected to be an even function of x. Using this fact, the odd 
functions of x of solution eq. 168 can be eliminated to yield 

. . R(a„nJ-a,X) 

W°(x) = F cosh(ax) cos(/?x) + G sinh(ax) sin(/?x) + — — 75 • l lby J 

A 11 A 2 2 A -(2 


The unknown A, are combined to form unknown constants F and G. By substituting the sol- 
ution form eq. 169 into the boundary conditions, constants F and G can be solved for. With F 
and G known, the displacement w?(x) is known. (Obviously with evenness of the solution be- 
ing enforced, the boundary conditions are enforced at one end of the cylinder only.) 

Numerical results for the three laminates of interest are illustrated in figs. 2-4. The cylin- 
der geometry and material properties used to obtain numerical results, in these figures and 
throughout, are as follows: The cylinder radius R is 10 in., the cylinder length is 30 in. 
(L/R = 3 , a short to intermediate length cylinder), and the cylinder thickness H is 0.080 in. 
(R/H = 125, a thin laminate). The material properties are 

E, = 20 Msi; E 2 = 1.3Msi; G 12 = 1.03Msi 

v 12 =--0.3; h (lamina thickness) = 0.005 in. (170) 

ai = -0.167 x 10“ 6 /°F; a 2 = 15.6 x 10“®/°F . 

In terms of laminate properties, for the (90 t /0 a ) T laminate 
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(171a) 


An = A 22 = 857,000;. A 12 = 31,400 Ib/in. 

Bn = 15,100 Ib/in.; = 457 Ib/in. 2 
N I = Nj = -259 Ib/in. /°F; M* = 3.94 lb/in. 2 /°F , 


For the (0 8 /90,) T laminate, 


A n = A 22 = 857,000; A 12 = 31,400 Ib/in. 

Bn = -15,100 Ib/in.; Dn = 457 Ib/in. 2 
= N J = -259 Ib/in. /°F; M* = -3.94 lb/in. 2 /°F 


For the (0/90)* laminate, 


A ii - A 22 = 857,000; A 12 = 31,400 Ib/in. 
Bn = 0; Dn = 532 Ib/in. 2 
Nx = Nj = -259 lb/in./°F; M r x = 0 . 


(171b) 


(171c) 


To neglect thermally-induced preloading effects, NJ , NJ and MJ are set to zero. Note that for 
all three laminates the elements of the A matrix are identical. The quantity B t1 is of opposite 
sign for the two unsymmetric laminates. Of course B„ = 0 for the symmetric laminate. Since 
the 0° layers are on the extreme outside and inside of the cylinder wall, the value of D,, is 
greater for the symmetric laminate. The quantities NJ and NJ are the same for all three lam- 
inates, whereas MJ is opposite in sign for the two unsymmetric laminates and zero for the 
symmetric case. The rearrangement of the 0° and 90° layers in the three laminates only in- 
fluence the out-of-plane, or bending, properties of the laminate. A temperature change of 
AT = — 280°F is used to represent the temperature change from a 350°F consolidation tem- 
perature to an ambient temperature of 70°F. 


The radial displacement, w?(x), as a function of length along the cylinder for the (90 8 /0 8 ) T 
laminate is shown in fig. 2. The displacements have been normalized by the laminate thick- 
ness while the axial coordinate has been normalized by the cylinder length. Several inter- 
esting characteristics of the thermally-induced response are illustrated in the figure. These 
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will be particularly interesting when compared with results for the (0 8 /90*) T laminate, to be 
discussed shortly. From the figure it is seen that over the majority of the length of the cylinder 
the displacement is uniform and inward, i.e., w < 0. Near the end, x/L > 0.4, there is a rapid 
change in the radial displacement, the displacement becoming more inward at the end of the 
cylinder (x/L = 0.5). This behavior is due to the unsymmetric nature of the laminate, in par- 
ticular, the thermally-induced moment generated within the laminate. Away from the ends of 
the cylinder the moment has no effect on the displacement. Thermal expansion in the x and 
6 directions, as well as the elastic properties in these directions, control the displacements 
away from the ends. At the ends of the cylinder there is no material to react the thermally- 
induced moment. The cylinder responds as if the ends were being subjected to an applied 
moment. The ends of the cylinder 'curl', producing a boundary layer effect. 

If the stacking arrangement is reversed and a (0,/90,) T laminate is considered, the cylinder 
deforms as shown in fig. 3. Interestingly enough, the central portion of the cylinder, away from 
the ends, moves radially inward the exact same amount as the (90,/0,) T cylinder. Near the 
ends there is the rapid change in radial displacement. In this case, however, at the end the 
cylinder deforms outward. This is opposite the situation for the (90,/0,) T laminate. The direc- 
tion the ends of the cylinder 'curl' is a function of the sign of the thermally-induced moment. 
The sign of the thermally-induced moment is a function of the stacking arrangement. The sign 
of the thermally-induced moment is opposite in the two cases, resulting in opposite bending 
deformations. The point to be made at this time is that before any loading is applied to the 
cylinder, there is a nonuniform radial displacement along the length of the cylinder. The lack 
of a uniform displacement, particularly since it is concentrated in a boundary layer near the 
ends, where the loads are applied, and particularly since the specific displacement is 
laminate-dependent and sign-sensitive, is an important point. 

For completeness, the thermal deformations of a (0/90) 4S laminate are illustrated in fig. 4. 
For this symmetric laminate the radial deformation due to the temperature change is uniform 
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along the length. The inward radial displacement is identical to the inward radial in the central 
portion of the (90,/0,) T and (0 8 /90,) T laminates. 

A note of interest: The deformation in the central portion of the cylinders, for any lami- 
nation sequence, is given by the particular solution, w^„(x). The response near the ends is 
controlled by the homogeneous solution. For the symmetric (0/90) 4S laminate the homogene- 
ous solution is zero. For the other two cases the homogeneous solution is a major influence. 
Since the particular solution is the same for each of the three laminates, the inward radial 
displacement is the same away from the ends for each of the three laminates. 

Of interest is the response that results when an axial compressive load is applied to the 
(90,/0 ,) t and (0 8 /90) T thermally-deformed cylinders. Is inclusion or exclusion of this thermally- 
induced preloading deformation important for predicting the load-induced response of the 
cylinder? These issues are discussed next. With or without the thermal preloading condition, 
calculating the response to an applied load is a more complicated problem. 

SOLUTION OF EQUATIONS FOR THE CASE OF THERMAL EFFECTS AND 

AN AXIAL END LOAD 

The combined problem of thermal and mechanical loads is governed by eq. 153 and 
boundary conditions eqs. 154, 155, or 156. As with the thermal-only problem, the solution 
consists of a homogeneous part and a particular part, i.e., 

w °(x) = w° omo (x) + wj art (x) . (172) 


The particular solution is again determined by inspection to be 


w part( x ) ““ 


^(AiiNfl — A 12 (N x t + N)) 


A 1i A 22 — A 


= w, 


part 


12 


(173) 
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Note that this particular solution is also independent of x. The homogeneous solution to 
this problem is taken to be of the form given by eq. 164. For this case, the characteristic 
equation is given by 




2A 12 B 11 - A 1Jt RN 


y+(-^ 


A 2 2 ~ A 12 


= 0. 


(174) 


Note the applied axial load appears in the characteristic equation. The level of applied load 
will be specified and hence all coefficients of k in the above equation are known. The four 
values of k are thus known at this point. 

The roots to the characteristic equation are given by 


^ 1 , 2 , 3, 4 — 


( A 11 RN-2A 12 B 11 ) ± 7(A^RN-2A 12 B 11 ) 2 -4(A 11 A 22 -A^)(D 11 A 11 -B^ 1 ) 

2(D 11 A 11 -B^)R 


(175) 

Though it is not totally obvious, there is an interesting character to the roots to eq. 175. This 
is due to the dependence of the roots on the level of applied load, N. The character of the 
roots can be examined by studying k * instead of k, i.e., 

2 (A^RN - 2A 12 B 11 ) ± ^(A^RN -2A 12 B-ii) 2 — 4(A 11 A 22 — A^ 2 )(P 11 A 11 — B^) 

71 : 2(D 11 A 11 - BiOR 


(176) 

The first important character to observe is that at a certain level of applied load the 
discriminant will be zero. In particular, the discriminant is zero when N = N*, N* being given 
by 
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(177) 



For a symmetric laminate B 1t — 0 and the determinant only vanishes for negative N*, a 
compressive load. For unsymmetric laminates, though it is unlikely, the determinant could 
vanish for positive N*. depending on the magnitude of A„ , A„ , , B„ , and D„. Here the 

negative sign will be assumed in eq. 177 and interest will focus on negative N*. 

Returning to eq. 175, for N = 0, the roots to A* are, in general, complex. The roots to A are 
thus of the form 

^i,2, 3 , 4= ± “ ±ip- (178) 

If the laminate is symmetric, N = 0 leads to two purely imaginary roots for A*. The roots for A 
are still of the form of eq. 178. 

For N increasing compressive the roots to A* remain complex, the roots for A being given 
by eq. 178. When N = N*. A 2 has negative real repeating roots. In particular, 

,i 2 

The four roots for A are then two repeating pure imaginary roots of the form 

^1,2, 3, 4= ± i/?. ± ip- (180) 

For N increasing in compression beyond N\ the roots for A * are distinct and negative real. 
The roots for A are of the form 

^ 1 , 2 , 3 . 4 = ± ifi it ± ip 2 - (181) 
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Because of these three different forms for the roots, depending on the value of N relative 
to N*. the functional form of the x dependence of the homogeneous solution depends on the 
value of N. This fact leads to a dependence on N of the deformed shape along the length of 
the cylinder. The amplitude of the deformed shape also depends on N. For the linear prob- 
lem, shape of the deformation does not depend on N. Only the amplitude of the deformation 
depends on N. This difference in the dependence on N between the linear and nonlinear 
cases is important. 

The functional dependence of homogeneous solution on x is as follows: 

For |N| < |N*|, from eq. 178 

»C,oM - A,8 < "' W< + A 2 e ( ’- | ' nx + A 3 e ( -“ +l » x + A 4 e ( -“- | » x . (182) 

This is the same form as for the thermally-induced preloading deformations. 

For |N| = | N* | , from eq. 180 

w°o,ro(x) = (A, + A 2 x)e i/?X 4- ( A 3 + A 4 x)e" ,/?X . (183) 

For |N| > |N*|, from eq. 181, 

WhomoM = A i cos^x) + A 2 sin^x) + A 3 cos (/? 2 x) (184) 

+ A 4 sin(/? 2 x) . 

Combining these solutions with the particular solution, eq. 173, and considering only the 
portion of the solution that is symmetric in x, the three forms are: 

For |N| < |N*|, 
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(185) 


w °(x) = F cosh(ax) cos(/?x) + G sinh(ax) sin(/?x) 
R(AiiNJ — A 12 (N x T + n)) 


+ ■ 


^ 11^22 — A 12 


For |N| = |N*|, 


w°(x) = F cos (fix) + G x sin(/?x) + 


^( A nNg — A 12 (N x + N)) 

a iiA 22 — A ? 2 


(186) 


For |N| > |N*|, 


w°(x) = F cos(/?.,x) + G cos(/? 2 x) + 


^(AiiNfl — A 12 (n x + N)) 

a iiA 22 — A ? 2 


(187) 


The constants F and G can be determined by the application of boundary conditions. The case 
of thermally-induced preloading deformations (N = 0), eq. 179, is a subcase of the solution 
given by eq. 185. 


With the above solution, the response of a cylinder to an axial end load in the presence 
of the thermally-induced preloading deformation can be determined as a function of N*. With 
the thermally-induced preloading the boundary conditions must be carefully interpreted. In 
particular, since the ends of the cylinder are already deformed due to the thermal preloading, 
the simple support conditions become 


W °( 2 ) — W °( ~2* ) 



dV 

dx 2 


( 

V R An 



In 

A ii 


(N -I- N x ) + M x 


= 0. 


(188) 


For a clamped boundary the conditions are 
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( 189 ) 


dw°(T) _ ^t(t) 

dx dx 


Any experimental set-up that is simulating simple support or clamped conditions would 
resist any end displacement beyond to the preloading value. Clamped conditions would pre- 
vent any rotation beyond the preloading value. Simple support and clamped conditions for the 
(90,/0 ,) t cylinder are depicted n figs. 5a and b, respectively. The situation for the (0/90) 4S cyl- 
inder is shown in figs. 6a and b. For unsymmetric laminates the load N* is applied eccen- 
trically relative to the midsurface of the cylinder wall. This could have a large influence on the 
response of the cylinder near the ends. 

NUMERICAL RESULTS FOR THE CASE OF THERMAL EFFECTS AND 
A COMPRESSIVE AXIAL END LOAD 

The primary reason for the current study is summarized in the figures to follow. Results 
from the closed-form solution of the radial displacement as a function of length along the cyl- 
inder for increasing compressive axial load levels are presented. The load levels examined 
are normalized by N*. Note that N* is independent of boundary conditions and thermal pre- 
loading. The value of N equal to N* has been interpreted as the value of N to cause buckling. 
Flugge [4] discusses this in connection with isotropic cylinders but the point is not pursued 
here. As will be seen, the response of the cylinder when N = N* suggests that N* may be 
associated with an instability. 

Simply-supported boundary conditions 

The deformed shape of a simply-supported (90,/0,) T cylinder for N = 0 is shown in fig. 7. 
As before, in this and subsequent figures the radial displacement has been normalized by the 
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thickness of the cylinder wall, H . The axial coordinate is normalized by cylinder length, L . 
Figure 7 is identical to fig. 2 except the radial displacement at the end of the cylinder is defined 
to be the location of zero displacement. In reality the radial displacement of the end of the 
cylinder is w T (l/2). This shift of coordinates is strictly cosmetic. If thermally-induced pre- 
loading deformations are neglected, the radial displacement is represented by the line 
w/H — 0. That the cylinder is assumed to be simply supported is irrelevant for this figure. 
The figure would be identical if the supports were considered clamped. Experimentally the 
cylinder would be mounted in the fixturing so as to effect the simple supports or the clamped 
supports. Since, in theory, no loads are involved in mounting the cylinder in the fixture, the 
preloading shape of the cylinder would not be influenced by the support conditions. Obvi- 
ously, when loads are applied to the cylinder, the support conditions become very important. 


The radial response of the simply-supported cylinder when the load is 10% of N* is shown 
in fig. 8. The response when the thermally-induced preloading deformations are incorrectly 
neglected in the analysis is indicated by the dashed line. This situation will be discussed. 
For the (90,/0,) T 


N* = —2460 Ib/in. ( 190 ) 

At the 10% level, the response due to the applied load is not too different than the thermally- 
induced preloading response, fig. 7. When the compressive load reaches 50% level, fig. 9, the 
applied load causes considerable response (note the vertical scale compared to fig. 8). The 
central portion of the cylinder has moved outward some, but the major effect is at the end of 
the cylinder. At the end of the cylinder the magnitude of the radial deformations has increased 
considerably, and though not detectable from the figure, the length of the boundary region has 
increased somewhat. The rapid change in magnitude at the cylinder end could clearly gen- 
erate significant bending strain. When the load is increased to 90% of N*. fig. 10, the boundary 
layer has clearly lengthened. The oscillatory nature of the deformations encompass one-half 
the cylinder. In addition, the magnitude of the deformations has become severe. In particular, 
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the cylinder well has displaced one wall thickness. Since there are high localized curvatures 
associated with this level of deflection. The level should cause concern. The oscillatory na- 
ture of the displacements is such that at points (e.g., x/L = 0.44) the cylinder radius is less 
than the undeformed radius, despite the fact that axial compression causes the cylinder, in 
general, to become larger in diameter. The stresses and strain that accompany these dis- 
placements are not reported cn here. 

As mentioned above, the dashed lines in figs. 8-10 represent the predicted response if the 
thermally-induced preloading deformations are not included in the analysis. As can be seen 
by examining the figures, neglecting the thermally-induced deformations results in a predicted 
displacement that is less than is predicted when the thermally-induced displacements are in- 
cluded. At N = 0.5N* the predicted peak displacements at x/L = 0.47 are different by a factor 
of two. At N = 0.9N* the peak displacements computed by neglecting the preloading thermal 
effects are about 40% lower than if preloading effects are included. 

The response of the simply supported (90,/0„) T cylinder when the applied load equals N' 
is shown in fig. 11. It is clear that the magnitude of the deformation is essentially unbounded 
at this load level. The deformations are oscillatory, the spatial frequency being given by eq. 
179. As stated previously, the present work will not investigate the relationship of this load 
level with the buckling load level. With the load increasing to 1.1N* the response of fig. 12 
results. The amplitude of the deformations has decreased considerably relative to the 
N = N* level. There are two basic frequencies to the oscillatory character of the response, 
these two frequencies being given by eq. 181. 

That the boundary layer increases in length as the load is increased is an important non- 
linear effect. For a linear analysis, the length of the boundary layer does not depend on the 
level of applied load. For this nonlinear problem, the length of the boundary layer and the 
level of applied load are related through the characteristic roots k, of eq. 179. In particular, the 
length of the boundary layer is related to the recriprocal of the real part of k, namely 1/a. The 
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values of 1 fa. as a function of the load level for the (90 8 /0 8 ) s laminate of figs. 7-10, and for the 
other laminates to be discussed shortly, are presented in Table 1. 


Table 1 

Boundary Layer Length, 1/a, for the Three Laminates 


N 

(90 8 /0 8 ) t 

(0 8 /90 8 ) t 

(0/90), 

0 

0.56 1/in. 

0.54 

0.71 

0.1N* 

0.59 

0.56 

0.74 

0.5N* 

0.79 

0.76 

1.00 

0.9N* 

1.77 

1.70 

2.23 


As can be seen, for the (90 8 /0 8 ) T laminate, the boundary layer triples in length as the load level 
changes from zero to 0.9N*. 

As a contrast to the response of the (90 8 /0 8 ) T cylinder, consider the response of a simply- 
supported (0,/90 ,) t cylinder. The characteristics for this cylinder are illustrated in figs. 13-18. 
The radial displacement at the condition of no axial load is shown in fig. 13, the displacement 
at the end of the cylinder again being taken as the zero displacement location. Figures 13 and 
3 differ only by a shift of origin. For the (0,/90,) T cylinder 

N* = -2681 Ib/in. (191) 

The sign of B„ is responsible for the difference between N* for the (90 8 /0,) T laminate and N* for 
the (0 8 /90 8 ) T laminate. 

As the load is applied to the (0 8 /90 8 ) T cylinder, the central portion of the cylinder moves 
outward, as it did with the (90 8 /0 8 ) T laminate. However, the boundary layer region develops 
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large inward deflections. At the location x/L - 0.47, the cylinder wall has moved in about one 
wall thickness when the load level reaches 90% of N*. Again the boundary layer region has 
grown considerably and developed rapid oscillations which could lead to large stresses. The 
length of the boundary layer, given in Table 1, is similar to that of the (90 8 /0 8 ) T laminate. It is 
interesting to compare figs. 16 and 10. In the boundary region, the rapid changes in the dis- 
placement of the (0 8 /90 8 ) t laminate are exactly opposite in sign to the rapid changes in the 
displacement of the (90./0,) laminate. The sign of B 11t which is opposite for the two laminates, 
is directly responsible for this feature. The response at N = N* is illustrated in fig. 17. The 
response is essentially unbounded at this load level. Comparing the response of this cylinder 
at N = N* with the response of the (90 8 /0,) T cylinder at N = N* , fig. 11, it is clear that the re- 
sponses of the two cylinders are exactly out of phase. This again is due to the sign of B 1V 
Since the values of N* are not the same for these two cases, the load is not the same in these 
two cases. The response for the case of N - 1.1N* is shown in fig. 18. The amplitude of the 
deformations has decreased relative to the N = N* level, and two basic frequencies are again 
evident. 

Since the unsymmetric laminates exhibit considerable radial deformation when com- 
pressed axially, it is of interest to study the radial deformations of a symmetric laminate. The 
response of the (0/90) 4S laminate is shown in figs. 19-24. The preloading radial deformations 
due to thermal effects are simply a spatially uniform displacement inward, as was shown in 
fig. 4. This uniform radial displacement is shown again in fig. 19 with the reference for zero 
displacement shifted, as has been done with the other two laminates. For this symmetric 
laminate the shifted preloading radial displacement coincides with the w = 0 line for all x. 
For the (0/90) 4S laminate 

N* = — 4270 Ib/in. ( 192 ) 

This is significantly higher than the values of N* for the other two laminates. 
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The radial displacement for a load of N equal 10% of N* is illustrated in fig. 20. The de- 
flection is everywhere outward. For the 50% and 90% levels, figs. 21 and 22, respectively, the 
deflections increase, and the boundary layer increases in length. The deflections at N* and 
1.1N* are shown in figs. 23 and 24, respectively. The deflections, though increasing with in- 
creasing N, do not reach the magnitudes that were observed for the two unsymmetric lami- 
nates. In addition, by comparing the values of N* for the three laminates, the loading on the 
symmetric laminate at 90% of N* is considerably higher than the loading on the unsymmetric 
laminate at 90%. Yet the displacements are a factor of 20 less with the symmetric laminate. 
Referring to Table 1, the boundary layer for this symmetric laminate is greater than for either 
unsymmetric laminate. When the thermally-induced preloading deformations are ignored in 
the analysis for the symmetric laminate, the response is predicted to be the same as when 
they are included. The lack of difference between the two situations can be attributed to the 
fact that the preloading deformations are spatially uniform for the symmetric laminate. There 
is no initial eccentricity at the ends of the cylinder where the load is applied. 

Clamped boundary conditions 

Since the simply supported unsymmetric laminates exhibited large, rapidly changing dis- 
placements near the end of the cylinder as the load level increased, it is legitimate to ask if 
similar behavior would be observed for other boundary conditions. It would seem that 
clamping the boundary would restrain the deformations considerably. The response of the 
(90,/0 ,) t laminate with clamped ends is illustrated in figs. 25-30. The case of N = 0 is shown 
in fig. 25, this being a repeat of figs. 2 and 7. This figure is repeated here for convenience. 
When a load of N equal 10% of N* is applied, fig. 26, the central portion of the cylinder moves 
outward and the deformations near the boundary begin to grow. As the load level is increased 
to 50% and 90% of N' , figs. 27 and 28, respectively, the boundary layer grows in length and 
the amplitude of the deformations in the boundary layer increase significantly, as they did for 
the simply-supported conditions. However, compared to the simply-supported case, figs. 8-12, 
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the clamped conditions sharply restrict the boundary layer deformations. At N = 0.9N , the 
deformations in the boundary layer for the clamped condition are one-fifth the deformations 
in the boundary layer for the simply-supported condition. The lengths of the boundary layers 
for the two support conditions are the same. The lengths of the boundary layers, given in 
Table 1, do not depend on the boundary conditions. At N = N*. the magnitude of the re- 
sponse of the clamped cylinder, fig. 29, is much less than for the simply supported case. If 
N* is associated with buckling, based on a comparison of figs. 11 and 29, it would have to be 
more closely associated with buckling for the simply-supported conditions than for the 
clamped conditions. For the clamped cylinder, at N = 1.1N*. fig. 30, the amplitude of the de- 
formation is greater than for N = N*. This was not the case for the simply supported cylin- 
ders, the response decreasing rapidly for N > N*. 

The dashed line in figs. 25-30 represent the predicted deformations when the thermally- 
induced preloading deformations are neglected for the clamped case. Neglecting the pre- 
loading deformations for the clamped case seriously under predicts the amplitude of the 
oscillating behavior in the boundary layer. At the 90% level the peak-to-peak amplitude of the 
oscillations of w/H is predicted to be over 0.1 if the preloading deformations are included, as 
they should be. If the preloading deformations are neglected, the peak-to-peak amplitude of 
the oscillations is quite small in comparison. 

Clamping the boundaries of the (0,/90 8 ) T laminate, figs. 31-36, produces a similar reduction 
in the amplitude of the oscillations in the boundary layer when compared to the simply- 
supported case. At the 90% load level the peak-to-peak amplitude of the oscillations is again 
about one-fifth the magnitude of the peak-to-peak amplitude of the oscillations for the simply- 
supported case. As with the comparison between the (90,/0,) T and the (0 8 /90 8 ) T laminates with 
simply supported boundaries, the sense of the boundary oscillations of the (0 8 /90 8 ) T laminate 
are opposite those of the (90 8 /0 8 ) laminate. At the x location where the radial displacement is 
most outward with the (90,/0 8 ) T , the radial displacement of the (0 8 /90 8 ) T is most inward. Also, 
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for the (0,/90 ,) t laminate, neglecting the preloading deformations leads to considerably differ- 
ent predictions for the response due to load. In fact, if the preloading deformations are neg- 
lected, the response of the (0,/90,) T and the (90,/0,) T are practically identical, as can be seen 
by comparing the dashed lines in figs. 25-30 with those in figs. 31-36. 

The effects of clamping the boundaries of the (0/90) 4S are shown in figs. 37-42, fig. 37 being 
trivial and a repeat of previous figures. As with the two unsymmetric laminates, clamping 
greatly reduces the amplitude of the oscillations in the boundary layer. At the 90% level and 
for simple supports, fig. 22, the peak-to-peak amplitude of the oscillations is about 0.05. 
Clamping, fig. 40, reduces the peak-to-peak level to less than 0.015. In addition, compared to 
the clamped unsymmetric laminates, the magnitude of the peak radial displacement of the 
clamped symmetric laminate is considerably less. Hence symmetric construction of the lam- 
inate reduces the displacements in the boundary layer independent of boundary conditions. 
Boundary conditions, however, have a significant influence on the displacements in the 
boundary layer. In this regard, it is important to note that boundary conditions have no in- 
fluence on the radial deformations in the central portion of the cylinder, away from the 
boundary layer. A close look at all the figures just discussed will show this to be the case. 
As mentioned before, the deformations in the central region are determined primarily by the 
particular solution to the problem. 

Before concluding this study, a comment should be made regarding the axial displace- 
ment. While attention has been focused on the radial displacement, the axial displacement 
can be easily determined. Using eq. 142 in eq. 149a and integrating, an expression for u°(x) 
can be obtained. By applying the condition that, for example u»(0) = 0, a complete expression 

for u°(x) results. The functional form of the integration depends, of course, on the value of N 
relative to N*. 
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CONCLUSIONS 

From the results presented, several important conclusions can be drawn. First, the 
inclusion of the thermally-induced preloading deformations in the analysis of response due to 
axial load has a significant inf uence on the predicted response. In general, if the thermally- 
induced preloading deformations are not included in the analysis, the predicted deformations 
are less. The specific amount less depends on the laminate, the boundary conditions, and the 
load level. A second conclusion is that unsymmetric lamination sequences result in much 
larger radial deformations than with symmetric laminates. The differences can be as much 
as an order of magnitude. Third, the sense of the lack of symmetry in the lamination sequence 
(i.e., the sign of the B SJ components) determines to a large degree the sign of the radial def- 
ormations. This is particularly true for the sign of the oscillatory deformation, near the 
boundary layer. Forth, for this nonlinear problem, the length of the boundary layer is a func- 
tion of the applied load level. In addition, the length of the boundary layer for the symmetric 
laminate is greater than the length of the boundary layer of either unsymmetric laminate. 
Finally, clamped boundary conditions suppress the radial deformation considerably relative to 
the simply supported case. The effects of asymmetry and thermally-induced preloading are, 
for the most part, independent of the support conditions. 
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Fig. 1 - Geometry and Nomenclature for Cylinder. 
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Fig. 3 - Room-Temperature Shape of the (0*/90,) T Cylinder (AT — 280°F). 
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Fig. 10 - Radial Deformations of the Simply-Supported (90 8 /0a) T Cylinder, N = 0.9N\ 


84 



ignoring 



including thermal deformations 



w/H 



0 0.1 0.2 0.3 0.4 

x/L 


Fig. 12 - Radial Deformations of the Simply-Supported (9OJ0Jj Cylinder, N 
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Fig. 16 - Radial Deformations of the Simply-Supported (0 8 /90 8 ) T Cylinder, N » 0.9N\ 
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Fig. 18 - Radial Deformations of the Simply-Supported (0,/90,) T Cylinder, N » 1.1 N*. 
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Fig. 22 - Radial Deformations of the Simply-Supported (0/90)* Cylinder, N ■ 0.9N\ 


96 








ignoring thermal deformations 


w/H 



101 


w/H 

0.2 



102 





including thermal deformations 




Fig. 30 - Radial Deformations of the Clamped (90,/0,) T Cylinder, N ■ 1.1 N*. 
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Fig. 31 - Radial Deformations of the Clamped (<W9 <W t Cylinder, N - 0. 
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Fig. 38 - Radial Deformations of the Clamped (0 
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Fig. 39 - Radial Deformations of the Clamped (0/90)^ Cylinder, N — 0.5N\ 
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Fig. 40 • Radial Deformations of the Clamped (0/90)4, Cylinder, N a 0.9N\ 
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